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HEAT TRANSFER TO BODIES TRAVELING AT HIGH SPEED IN THE UPPER ATMOSPHERE 

By Jackson R. Staldeb and David Jukoff 


SUMMARY 

A general method has been developed , using the methods of 
kinetic theory, whereby the surface temperatures of bodies can 
be calculated for steady flight at any speed in a rarefied gas. 
The particular solution was made for aflat plate; however, the 
calculations can be easily extended to bodies of arbitrary shape. 

It was found that the aerodynamic heating problem in the 
absence of solar radiation, that is, for the case of nocturnal 
flight, becomes of negligible importance at altitudes of 125 miles 
and higher and up to steady flight speeds of 36,000 feet per 
second. The effect of solar radiation, for the case of daytime 
flight, becomes increasingly important as the flight altitude is 
increased. At an altitude of 150 miles and higher, solar radia- 
tion is the predominating factor that determines skin tempera- 
ture. Owing to the strong effect of solar radiation on skin 
temperatures at high altitudes, the desirability of nocturnal 
flight is indicated in order to minimize skin temperatures. 

In order to maintain low skin temperatures, it was found 
that the angle of inclination of the body with respect to the 
flight path should be kept as small as possible. This may be 
accomplished in practice by designing the body to be finely 
tapered and by flying the body at small angles of attack. 

It is pointed out that skin temperatures may be reduced by 
insuring thermal contact between portions of the skin inclined 
at positive and negative angles with respect to the flight path. 
As much surface as possible should be inclined at negative 
angles. Practically, this may be accomplished by boattailing 
the body. 

In the event that an internal skin-cooling system is employed, 
it is shown that the rate of internal cooling must be of the same 
order of magnitude or greater than the rate at which heat is lost 
naturally by emitted radiation. If the cooling rate is below 
the natural radiation rate, cooling has little effect upon skin 
temperatures . 

It is shown that, in the case of a missile designed to fly over 
a wide range of altitudes and speeds, it is desirable to make the 
emissivity of the skin as high as possible. This conclusion, 
however, is based upon a skin surface for which the emissivity 
is independent of the wave length of the emitted and absorbed 
radiant energy. 

A possible method of reducing surface temperatures is indi- 
cated by the decrease in skin temperature which accompanies a 
decrease in thermal accommodation coefficient. This phenome- 
non may be used to advantage if it is possible to decrease the 
accommodation coefficient by altering the surface characteristics 
of the skin. 


INTRODUCTION 

Missile or aircraft flight at extremely high altitudes intro- 
duces many new problems to the designer. One of the more 
important problems that arises is due to the skin tempera- 
tures that are attained at very high velocities. The severity 
of the problem is commonly exemplified by consideration of 
the disintegration by combustion of most of the meteors 
that enter the earth’s atmosphere at high speed. 

In reference 1, Tsien has summarized the meager amount 
of research that has been done in the field of rarefied gas 
kinetics as applied to aeronautics. It is shown that the 
general field of aeronautics may be divided into three major 
regimes, and that the division lines, although poorly defined 
at present due to lack of experimental data, are marked by 
the ratio of the mean free molecular path length to a char- 
acteristic body dimension. 

The three main regimes of aerodynamics, or gas kinetics, 
may be described very briefly as: 

1. The region of conventional aerodynamics in which the 
mean free path length is negligible in comparison with the 
body size. 

2. The region of free-molecule phenomena in which the 
mean free path is large with respect to the body size. 

3. The so-called “slip flow” regime, intermediate between 
the conventional aerodynamics and free-molecule regimes. 

In conventional aerodynamic theory, the first basic as- 
sumption that is made is that the fluid under consideration 
may be treated as a homogenous continuous medium. This 
assumption certainly should become invalid at very high 
altitudes owing to the extreme rarefaction of the atmosphere. 
A measure of the atmospheric density is furnished by the 
average length of path (mean free molecular path) that is 
traveled by a molecule between successive collisions. The 
mean free molecular path has been computed by Johnson and 
Slack in reference 2 for altitudes up to 150 miles. The result 
of their computations is shown in figure 1. The existence 
of a large mean free path should not be construed as an indi- 
cation that the number of molecules per unit volume is small, 
in the ordinary sense of the word. For example, at a reduced 
density at which the mean free path is 10 feet, there are still 
10 13 molecules contained in a cubic inch. 

The work of this report will be confined to a study of the 
heat-transfer processes occurring in the region of free 
molecule flow. The lower limit of the ratio of mean free 
path to a characteristic body dimension at which the derived 
theory becomes invalid cannot be stated at this time for 
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lack of experimental evidence. The numerical computations, 
however, will be confined to altitudes of 75 miles and higher. 
The estimated mean free path at 75 miles altitude is 1 foot. 

The specific purpose of the present report was to develop a 
method whereby the surface temperatures of bodies travel- 
ing at any speed in a rarefied gas could be determined. It 
should be pointed out that, to the authors’ knowledge, no 
experimental work has been done in high-speed rarefied gas 
streams. Consequently, the numerical values obtained 
should be regarded with some reservation until experimental 
evidence is accumulated concerning the validity of the as- 
sumptions that are made in the analysis. 

SYMBOLS 

constant of integration in Maxwell’s equation 
total velocity, feet per second 

components of total velocity, feet per second 

base of natural logarithms, 2.7183 . . . 
energy of incident molecules, foot-pounds per 
square foot, second 

energy of molecules at plate temperature, foot- 
pounds per square foot, second 
energy of re-emitted molecules, foot-pounds per 
square foot, second 

energy of rotation, foot-pounds per molecule 
translational energy of molecules striking front of 
plate, foot-pounds per square foot, second 
translational energy of molecules striking rear of 
plate, foot-pounds per square foot, second 
vibrational energy, foot-pounds per molecule 

Maxwell’s distribution function, dimensionless 

gravitational acceleration, assumed constant at 
32.2, feet per second squared 
number of degrees of freedom of motion 
normal solar radiation at top of atmosphere, 93.4 
foot-pounds per square foot, second 
Boltzmann constant, 5.66X10 -24 foot-pounds per 
degree F per molecule 
mass of one molecule, slugs per molecule 
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C x 

Cy 

c 2 

e 

E t 

E P 

E r 

E R 

E, 

E/ 

E„ 

J(V) 1 

r ( c ) j 
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J . 
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No 
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T t 

T P 

Tr 

u 

U x 

Uy 

Uz 

V 

v m 

v z 

Vy 

Vz 
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i6 
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$(a) 

x 

x' 

* 
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molecular weight, pounds per pound-mole 
number of molecules striking front of plate, mole- 
cules per square foot per second 
number of molecules striking rear of plate, mole- 
cules per square foot per second 
number of molecules per unit volume of gas 
Avogadro’s number, 2.73 X10 26 molecules per 
pound-mole 

heat removed from skin by an internal cooling 
system, foot-pounds per square foot, second 
Universal gas constant, 1544 foot-pounds per 
pound-mole, °F absolute 

incident radiant energy, foot-pounds per square 
foot, second 

re-emitted radiant energy, foot-pounds per square 
foot, second 

molecular speed ratio (jy - ) 

surface area, square feet 
absolute temperature, °F absolute 
temperature of incident molecules, °F absolute 
temperature of body skin, °F absolute 
temperature of re-emitted molecules, °F absolute 
mass velocity, feet per second 

components of mass velocity, feet per second 

thermal velocity, feet per second 

most probable molecular speed, feet per second 

> components of thermal velocity, feet per second 


space coordinates 

accommodation coefficient, dimensionless 
reciprocal of most probable molecular speed, 
seconds per foot 
emissivity, dimensionless 
angle of incidence, degrees 

Stefan-Boltzmann constant, 3.74 X10~ 10 foot- 
pounds per square foot, second, °F 4 absolute 


2 r 

error function, — = e 

VvJo 


dx, 


dimensionless 


dimensionless group defined by equation (7) 
dimensionless group defined by equation (9) 
dimensionless group defined by equation (14) 
dimensionless group defined by equation (19) 


ANALYSIS 

In order to determine the skin temperature of a body mov- 
ing at constant speed through a rarefied gas, it is necessary, 
fundamentally, to make an energy balance on a section of 
the skin. Energy may be added to or subtracted from the 
skin by three distinct processes: 

(а) Molecular energy transport to and from the skin. 

(б) Radiant energy transport to and from the skin. 

(c) Energy added to or removed from the skin by processes 
occurring within the body. 











. £^3Sfe' 


Av* 


HEAT TRANSFER TO BODIES TRAVELING AT HIGH SPEED IN THE UPPER ATMOSPHERE 




Thus, we may write a simple energy balance on a section of 
the skin as 


y 



Ei+R i =E T -\-R T J rQ ( 1 ) 

Each term in this equation can be considered separately. 

MOLECULAR ENERGY TRANSPORT TO THE BODY 

It is a matter of common knowledge that any body sur- 
rounded by a gas is constantly bombarded by molecules of 
the gasi Tt'is tins molecular bombardment, in fact, that 
causes a pressure to be exerted on a body or on the walls of 
a container in which the gas is confined. It is evident that 
the molecules possess, by virture of their mass and velocity, 
a certain kinetic energy. Thus, a body immersed in a gas is 
constantly receiving and emitting energy from the bombard- 
ment of impinging molecules and the re-emission of the same 
molecules from the body surface. If it were considered that 
the molecules carried no energy except that due to transla- 
tional velocity, then, we should, by examination of their 
incident and re-emitted velocities, be able to determine the 
net increase or decrease of the kinetic energy of the re-emitted 
gas and, hence, determine the molecular energy added to or 
subtracted from the body. In the ordinary case, however, 
the impacts of the molecules with the body become obscured 
by the collisions of the molecules among themselves and the 
problem becomes exceedingly complex if treated on a micro- 
scopic basis. On the other hand, if we consider a body 
immersed in a gas that is of low enough density so that 
collisions among molecules in the vicinity of the body are 
rare in comparison with collisons with the body, the problem 
becomes susceptible to analysis by the conventional methods 
of kinetic theory. 

If it is assumed that the gas through which the body is 
moving is so rarefied that the motion of the impinging 
molecules is essentially unaltered by collisions with reflected 
molecules, then it follows that the velocity distribution of the 
impinging molecules will be that of a gas in thermal equilib- 
rium, the Maxwell distribution. The Maxwell distribution 
law for a gas at rest states that of a total of N molecules per 
unit volume, the fraction that has velocity components lying 
between the values, V x to V x +dV x , V v to V v +dV,„ V z to 
V z -\-dV z , is given by 

f(V x V v V z )dVJV y dV x =Ae-^ v ^ v ‘ 2+v ^dV x dV,dV z (2) 



plate. All of the molecules that strike dS during a time 
interval dt that have velocities lying in the range C x to C x - f- 
dC x , C v to Cy+dCy, and C z to Cz+dCz will lie in a cylinder 
having dS as a base and with a length C dt in the direction of 
C, where C— VC J 2 +C'„ 2 +C'z 2 . The volume of the cylinder 
is C x dS dt. The number of molecules of this kind contained 
in the cylinder is then/' {C x , C v , Cz) NC X dS dt dC x dC v dC z 
where/' ( C x , C y , C z ) is Maxwell’s distribution function for the 
total velocities. Then, to find the total number of molecules 
striking dS per second, we must integrate j' (C x , C y , C z ) 
NC X dS dC x dC y dCz between the limits C x — 0 to °o, C v — 
— co to <», and Cz ~ — 00 to ®> . Molecules having a nega- 
tive velocity in the x direction cannot strike the front side 
of the plate. Then, 

n=ArJ” J” J fl ° (C x , Cy, Cz) C x dC x dC y dC, (3) 


»• 

i 




I 



where V x , V v , V z are the random thermal velocities with 
respect to a coordinate system x, y, z fixed in space. v. 

To determine the energy transported to the body by the 
impinging molecules, it is necessaiy first to determine the 
number of molecules that strike the body. This is done as 
follows: 

Assume that the body is stationary and that the gas has a 
macroscopic velocity TJ with components U x U„, U z with 
respect to the fixed coordinates. Assume the body to be a 
flat plate and fix the plate at the origin of x, y, z, so that the x 
axis is normal to the plate. The mass velocity U makes an 
angle 0 with the plate. Figure 2 is a diagram of the co- 
ordinate system. 

The total velocities may be written as the sum of the ther- 
mal and mass velocities, thus, C x — U x +V z , C y —U y +V v , 
and C z =U z -\-V z . Consider an elemental area dS of the 


The distribution function for the thermal velocities with a 
superimposed mass velocity may be written as (cf., reference 
3, p. 45) 

/' ( C x , C v , C x ) =Ae~P 1 


Thus for a unit area, 


n=NA rss \-mc,-v^+,Cz-u^c.-u.)» CxdCxdCvdC! 

(4) 

This equation can be integrated; making use of the relation 


B 3 

A—- 1 ^ the result is 

7T 3 ' 2 


( TJ ) 

Hlvs+f [I +* ( 0 p ->| ® 
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The factor /3 is related to the most probable molecu lar speed 
V m as fi=l/V m , while V m is given by V m =^2 gRT/M. Also 
U x may be written as TJ sin 6. Equation (5) may then be re- 
written as 




i(, + -yjtc s sin 8 [l+$(s sin 0)]} 


( 6 ) 


This is the result obtained by Sanger in reference 4 and re- 
derived by Tsien in reference 1. A plot of the dimension- 
less group x, defined as 

X={e~ $2 aia!! e + Vir S sin 8 [l + $(s sin 0)]} (7) 

is shown in figure 3. 


16 
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0 .8 1.6 2.4 3.2 4.0 

a sin 8 

Figure 3.— The dimensionless quantity x . as a (unction o( a sin 6. 
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To determine the number of molecules striking the rear 
side of the plate n' , equation (4) is again integrated except 
that the limits of G x are changed to C x =— m to 0. The 
result is 


The group 


_ NVrn 

~2 y+ 


l'={e 


| e - S 2ain2(,_ ^ g [l— $(s gi n #)]} (8) 


■ -\pK s sin 6 [1— <t>(s sin 0)]} (9) 


The total kinetic energy of the molecules striking the front 
side of the plate is obtained by integrating equation (10) 
over the limits C x —0 to Cy— — “> to ”, and C z = — 00 
to co . Then, for a unit area, per unit time, 


is plotted in figure 4. 

To determine the kinetic energy incident upon the plate 
due to translation, Sanger in reference 4, and Tsien in refer- 
ence 1, multiply the mass of gas striking the plate nm by 
the sum of the kinetic energy per unit mass due to the mass 
motion of the gas U 2 I2 and the average kinetic energy due to 
thermal translation which may be shown, for a static gas, 

However, C 2 =C x 2 +C y 2 +C z 2 , so that equation (11) may be written as 

1 Jo C x {C x 2 +C„ 2 +C 2 )e-^ w.-o^Hci t -u.wric.-u.n dC x d,C v dC, 

1 



0 .4 .8 1.2 1.6 2.0 2.4 

s sin 8 

Figure 4.— The dimensionless quantity as a function of s sin 0. 


3 R T 

to have the value ^his procedure, however, is an 

oversimplification of the problem as introduction of the term 
3 RT 

2 § M assumes that the kinetic energy due to thermal 

velocity may be added to the kinetic energy due to mass 
velocity. That this is not permissible is explained by the 
fact that the total kinetic energy depends upon the square 
of the total velocity, hence the components of the kinetic 
energy due to the velocity components of the total velocity 
may not be added linearly. Symbolically, this may be 
expressed in simple form as 

mV 2 mU 2 _ x m(U-\-V) 2 
~2 1 2 ^ 2 

A more exact calculation of the incident translational 
energy can be made as follows: 

The kinetic energy per molecule contained in a cylinder 
placed on dS is mC 2 !2. The kinetic energy of those molecules 
in the cylinder with velocity components between C x to 
C x +dC x , Cy to Cy+dCy, C z to C z +dC z is 
m 


NC 2 C x f ( C x ,C v ,C z )dS dt dC x dC v dC z 


( 10 ) 


E <=™ N J o ” C 2 CJ' (C X ,C V 


C 2 ) dC x dCy dC z (1 1) 


TO AT . 

E,= g NA 


( 12 ) 


The result of a tedious integration is 


E t =nm 





l + ^pU x [l+<Hl3U x )]e^ )_ 

where n is given by equation (6). Using the relations, 


(13) 
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Z WT 


0 = 


'V, 


> U x = U sin 


and defining a dimensionless group \p as 
*=1 


1+1 V ir s s i n 0 [l+4 , (s sin 0)] e s2sln!<! 

Jj 


1 + VV£ sin d [l+$(s sin 6)] e : 


,S 2sln2 0 


(14) 
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sy Equation?*(13) may be written as, 



E--n (15) 

When the mass velocity U is taken as zero, equation (15) 
reduces to 

E,~2nkT (16) 

which~is~the usual result obtained for the translatory kinetic 
energy crossing ah imaginary plane drawn in a static gas. 
When s sin 9 is large, the equation approaches 

E t =n(^P+\kT^ (17) 

A similar computation yields, for the incident translational 

energy on a unit area of the rear side of the plate, 

E t '=n' (18) 

where 

1 Vtt s sin 0[1— $(s sin fl)]e f>2sin2< ' 

^' = 1-1 — (19) 

1 — Vv s sin 0[1— $(s sin 0)]e s sin29 

The value of E,' , of course, reduces to 

E! = 2n'kT= 2nlc T (20) 

for U equal to zero and approaches E,' — 0 for large values of 
s sin 8. The dimensionless groups ^ and \ p' are shown in 
figures 5 and 6, plotted as functions of s sin 9. 



If the gas is considered to be composed of hard-sphere 
molecules, then, for a monatomic gas, the total kinetic energy 
incident upon the plate is given by equations (15) and (18) 
for the front and rear sides, respectively. However, if the 
gas is composed of diatomic dumbbell molecules, each mole- 
cule may carry an additional amount of energy due to rota- 
tion and vibration. By the principle of the equipartition 
of energy, the maximum amount of energy carried by a 

j 

molecule is 2 kT where j is the number of degrees of freedom 

of motion. A rigid dumbbell molecule may rotate about 
two mutually perpendicular axes, hence the maximum aver- 
age kinetic energy of rotation per molecule is 



E R =kT (21) 


It is conceivable that collisions between dumbbell molecules 
could cause the atoms composing the molecule to vibrate or 
oscillate along the line joining their centers. In this case, 
the energy of vibration is composed of kinetic and potential 
energies and the vibrational energy is equally distributed 
between the kinetic and potential forms. Thus, the total 
energy due to vibration is 

E v —2 (i JcT^—kT (22) 

where the factor 2 accounts for the distribution of energy 
between the kinetic and potential forms. 

Some question may arise as to the correct procedure to 
follow when considering a gas composed of two or more 
constituents. It can be shown that the results that have 
been obtained for a homogeneous gas, may be extended to 
mixtures of gases by calculating n, E,, E r , and E, separately 
for each constituent. This is a result of the fact that, in a 
mixture of gases in equilibrium, each gas has the same velocity 
distribution that it would have if the-other constituents were 
not present. 

The result obtained in equations (15), (16), (18), and (20) 
shows that, in effect, both front and rear sides of the plate 
selectively collide with a group of molecules having an average 
thermal velocity in excess of the average thermal velocity 
of the molecules of the gas considered as a whole. This is 
true whether the plate is moving or stationary. As the plate 


i 


y 
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velocity increases, the average thermal velocity of those 
molecules striking the front side approaches the average 
thermal velocity of the molecules of the gas considered as a 
whole, while the average thermal velocity of those molecules 
striking the rear side departs farther from the average thermal 
velocity of the molecules of the main mass of gas. Because 
of this selective interception, the plate is struck on both 
front and rear sides by molecules that possess higher-than- 
average translational energy associated with their higher- 
than-average thermal velocity. It is pertinent to inquire if 
this does not mean, also, that this nonaverage group of 
molecules carries nonaverage rotational and vibrational 
energies. It is pointed out in reference 5 that Maxwell’s 
distribution law allows the assumption to be made that 
rotational and translational energies are completely indepen- 
dent, so that a molecule with a high translational energy is 
just as likely to have a low as a high rotational energy. In 
the case of vibrational energies, it is shown that those mole- 
cules that have higher translational velocities also carry 
higher vibrational energies. However, for reasons which 
will presently be discussed, the change in the vibrational 
energy is neglected in the energy balance equation, hence we 
need not be concerned with its magnitude in this simplified 
analysis. 

MOLECULAR ENERGY TRANSPORT FROM THE BODY 

In order to determine the energy carried from the skin by 
molecules returning to the main mass of gas, it is necessary 
to consider the mechanism of energy transfer and molecular 
reflection at the surface of the body. A complete discussion 
is beyond the scope of this paper; the reader is referred to 
references 5 and 6 for a detailed description of the phenomena 
involved. Nevertheless, the more important concepts are 
outlined below because an understanding of them is necessary 
in succeeding computations. 

Two main types of reflection may occur when a molecule 
strikes a surface. The first and simpler type is specular 
reflection. In this case, a molecule leaves the surface with 
its normal velocity component reversed, retaining unchanged 
its tangential component. Any other type of rebound from 
a surface is known as diffuse reflection. It is possible to 
devise many subtypes of diffuse reflection; however, only 
a limited number are physically probable or even thermo- 
dynamically permissible. Insofar as direction of reflection 
or re-emission is concerned, it can be shown that the proba- 
bility of a molecule leaving the surface at a specified angle 
is proportional to the cosine of the angle with respect to the 
normal. All directions about the normal are equally prob- 
able, regardless of the direction of impingement. Several 
hypotheses may be advanced concerning the speed at which 
diffusely reflected molecules are emitted. The most physi- 
cally likely assumption is that the speeds are grouped in a 
Maxwellian distribution corresponding to a temperature 
which is not necessarily that of the surface. 

It appears logical that diffuse reflection would predominate 
as, in reality, an impinging molecule collides, not with a 
smooth hard surface, but with an individual molecule or 
atom of the body material. It would seem likely that col- 
lisions of this type would obliterate, insofar as directions are 


concerned, any prior directional history of motion, especially 
if more than one collision occurred at the surface. 

As yet, the temperature at which the molecules leave the 
surface has not been specified. In general, the molecules 
that leave the body will not possess an average velocity cor- 
responding to the plate temperature. In the case with which 
we are concerned, namely, that of a cooled plate, the molecules 
will leave at some temperature which is higher than the plate 
surface temperature. This inequality of temperature (or 
energy levels) is defined by introduction of the accommoda- 
tion coefficient, defined as 


The accommodation coefficient may be considered as a term 
that accounts for the inability of the impinging molecules to 
adjust themselves to the plate temperature during the time 
they are in contact with the plate. As pointed out by Loeb, 
in reference 5, it appears that the action of surface roughness 
is different for directional effects and energy or temperature 
effects. Apparently, one or two collisions with the molecules 
of the plate is sufficient to erase the directional history of a 
molecule; however, many collisions are required to attain 
temperature equilibrium with the plate. In the definition 
of the accommodation coefficient given by equation (23), the 
effects of specular and diffuse reflection on the energy ex- 
change are included in the value of a. 

The value of the accommodation coefficient has been 
measured for several static gases in contact with many dif- 
ferent materials. Values reported range from less than 0.1 
to 1.0, as shown in references 7 to 12. Reference 13 presents 
a method for correlating the effect of pressure upon accom- 
modation coefficients. Apparently, no measurements have 
been made of the accommodation coefficient in a high-speed 
gas stream. It would appear probable, however, that the 
accommodation coefficient would depend upon both the 
speed of the gas stream and the angle of the surface with 
respect to the stream. 

The accommodation coefficient, if defined rigorously, 
would have different values for translational, rotational, and 
vibrational energies. However, a series of experiments per- 
formed by Knudsen and described in reference 5 showed that 
the accommodation coefficients were numerically identical 
for rotational and translational energies. The vibrational 
components of the internal energy require much more time to 
adjust to new values than do the rotational and transla- 
tional components. A discussion of the time of adjustment 
of vibrational energy for gases may be found in reference 14. 
In succeeding calculations, it will be assumed that the vibra- 
tional component of the internal energy is unchanged during 
the short time interval that a molecule is in contact with the 
body. By the same line of reasoning, dissociation of poly- 
atomic molecules is neglected, since dissociation is a func- 
tion of the intensity of the molecular vibrations. These as- 
sumptions are serious in that vibrational energy and the 
energy absorbed during dissociation would, if present, have 
a marked effect on the plate temperature. If the length of 
time a molecule is adsorbed on the plate is equal to or 
greater than the time required for fully developed vibration 



CALCULATIONS 


to be excited, then dissociation would probably occur at the 
high plate temperatures calculated for high-speed flight. A 
more complete analysis awaits experimental data on ad- 
sorption times and accommodation coefficients in a flowing 
gas. 

To calculate the energy abstracted from the plate due to 
the return of molecules to the stream of gas, use is made of 
equation (23) rewritten as 

E r —Ei(l — a) + aE p (24) 

If E v is the energy of the molecules if they were at plate 
temperature, then equation (24) can be written for a mon- 
atomic gas that obeys the perfect gas laws as 

E T =E,(l — a)+2ankT P (25) 

and for a diatomic gas that obeys the perfect gas laws, 
neglecting vibrational energy changes, as 

E r = {E t + nk T t ) (1 ■ — a) + 3 ank T v (26) 

RADIANT ENERGY EXCHANGE WITH THE BODY 

Any body is constantly receiving and emitting radiant 
energy. In the steady state, a balance is soon reached where 
the emitted energy is equal to the incident energy; the body 
temperature adjusts itself to a value such that this equi- 
librium is maintained. 

In the present case, the body emits radiation at a rate 
given by 

Br=t*T* (27) 

The body may absorb radiation by several processes. 
Probably the most important is that due to direct solar 
radiation . The mean value of the solar radiation at the top 
of the earth’s atmosphere, the solar constant J 0 is shown in 
reference 15 to have a maximum value of 93.4 foot-pounds 
per second, per square foot for a plate normal to the incident 
radiation. The body may also absorb radiation from con- 
stituents of the atmosphere and from the earth’s surface. 
Most of the atmospheric gases do not emit appreciable radi- 
ation; however, water vapor, ozone and carbon dioxide emit 
strongly at certain wavelengths. In the interest of simplic- 
ity, radiation from the atmospheric constituents and from 
the earth’s surface will be neglected and, where incident 
radiation is considered, it will be assumed to be due solely 
to solar radiation, having its previously mentioned maximum 
value of 93.4 foot-pounds per second, per square foot. The 
amount of this energy that is absorbed is then 

R t =tJ 0 ( 28 ) 

provided that gray body absorption is assumed; that is, the 
emissivity is constant for all wavelengths and temperatures. 

ENERGY EXCHANGES OCCURRING WITHIN THE BODY 

It is conceivable that, in a missile designed for sustained 
flight at very high speeds, a skin-cooling system could be used 
to maintain structure temperatures within specified limits. 
Disregarding a discussion of the possible cooling systems, we 
may simply assume that the cooling system removes heat 
at a rate Q foot-pounds per second, per square foot. 

868900—50 2 


Having outlined the processes by which energy is ex- 
changed with the body, it is now possible to calculate tho 
temperature of a plate moving at any speed through a 
rarefied atmosphere. 

METHOD OF CALCULATION 

We may substitute equations (25) or (26), (27), and (28) 
in equation (1) and obtain, for a monatomic gas, where 
E t =Et 

aE,-2ankT J ,+eJ o -eaT P i -Q=0 (29) 

and for a diatomic gas, where E i =E t J r nE Il 

ct(E,-\-nkTi ) — 3ankT p -}-eJ„ — t<xT P i — Q = 0 (30) 

In equations (29) and (30) , the value of E t may be taken from 
equation (15) and the value of n from equation (6) for the 
front side of the plate. For the rear side of the plate, equa- 
tions (18) and (8) are substituted for values of E and n' , 
respectively. 

For a mixture of monatomic and diatomic gases, the inci- 
dent energies of the components of the gas are added together, 
and the re-emitted energies of the components of the gas are 
added together, so that the energy equation assumes the form 

T,E l -XE r +eJ 0 -€<rT 1> i -Q=Q (31) 

To determine the temperature of a plate with no cooling 
except that due to radiation, the factor Q in equations (30) 
and (31) is set equal to zero and the equations are solved for 
values of T p . 

To determine the effect of internal cooling on the plate 
temperature, arbitrary values of Q are substituted into 
equations (30) and (31) and the equations are again solved 
for values of T v . 

SCOPE OF CALCULATIONS 

Calculations were made to determine the temperature of a 
plate traveling at velocities of from zero to 36,000 feet per 
second at altitudes of from 75 to 190 miles. A speed of 
36,000 feet per second is approximately the minimum speed 
required to escape from the earth at these altitudes. The 
plate, in the first case, was assumed to be insulated between 
the front and rear surfaces so that the two surfaces could be 
treated independently. In the second case, the front and 
rear sides were assumed to be in perfect thermal contact ; that 
is, the temperature of the front side of the plate was equal to 
that of the rear side of the plate. In these two sets of calcu- 
lations, it was assumed that solar radiation was absent. 

The above calculations were then repeated, for one altitude, 
and angle of incidence (or body angle) with arbitrarily 
assumed values of internal cooling and with solar radiation 
absent. 

A fourth set of calculations was made for several altitudes, 
considering that the plate was uncooled and that solar 
radiation was present. 

In the first four sets of calculations the values of the 
accommodation coefficient and the emissivity were assumed 
to be. constant at unity. This assumption simplified the 
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computations and, at the same time, reasonably approxi- 
mated physical conditions so that the trends of skin temper- 
ature shown should be correct in order of magnitude. 

A fifth and sixth set of calculations were made to show the 
effect of emissivity and accommodation coefficient on skin 
temperature at several altitudes and flight speeds. 

The atmospheric properties used in the calculation were 
taken from references 2 and 16. Reference 16 was used to 
obtain the properties of the atmosphere at 75 miles; whereas 
reference 2 was used for the higher altitudes. There is some 
difference in the assumed atmospheric compositions between 
the tables of references 2 and 16. Warfield, in reference 16, 
assumed that the composition by weight at all altitudes was 
constant at its sea-level value, except for dissociation of 
oxygen. The amount of dissociation was assumed to in- 
crease with altitude, complete dissociation occurring above 
altitudes of 75 miles. On the other hand, Johnson and Slack 
assumed sea-level composition at all altitudes, neglecting 
dissociation entirely. In order to compare values of skin 
temperature at the several altitudes for which computations 
were made, it was arbitrarily decided to use the values of 
temperature and density given by Johnson and Slack and to 
assume complete dissociation of the assumed constant sea- 
level percentage of oxygen. The following table lists the 
values of atmospheric density and temperature that were 
used in the computations. 


Altitude 

(miles) 

Temperature 
(°F absolute) 

Density 

(lb/ft 3 ) 

75 

680 

2. 55X10-8 

100 

1025 

1. 05X10- 9 

125 

1385 

1.00X10- 10 

150 

1735 

1. 75X10-» 

190 

1990 

1.61X10" 14 


DISCUSSION 

INCIDENT TRANSLATIONAL ENERGY 

In the computation of the incident kinetic energy due to 
molecular translation, it was shown that both the number of 
molecules striking unit area of the body and the resultant 
kinetic energy incident upon the body were functions of 
the term s sin 9. An examination of the graphs of the 
dimensionless groups x and x'> ’A an d 'P' , which govern the 
number of molecules striking a section of the body and the 
magnitude of the incident translational energy, reveal 
several interesting features. With respect to the number of 
molecules that strike a body on the front side, it can be seen 
from figure 3 that this number increases practically linearly 
with flight speed for a given plate or body angle. On the 
other hand, as shown in figure 4, the group x' that deter- 
mines the number of molecules that strike the rear side of a 
body, decreases rapidly with flight speed and is practically 
zero for values of s sin 6 greater than 2. 

From equations (15) and (18), it may be seen that the 
total incident translatory kinetic energy term is composed 
of two parts: The first due to mass motion of the gas rela- 
tive to the body and the second due to the apparent trans- 
latory molecular motions. In the case of the front side of 
the body, the apparent molecular energy term can be seen, 
from figure 5, to increase rapidly at first with increasing 


values of s sin 9 and then to approach a constant value of 5/2 
at values of s sin 6 greater than 2. On the rear side of the 
same body, the molecular energy term continues to decrease 
rapidly as s sin 6 increases. 

A simple physical interpretation of the trends noted 
above can be given. First, it is necessary to understand 
that with a Maxwellian type of velocity distribution, most 
molecules have speeds close to the average molecular speed. 
A very small fraction has speeds much greater than the 
average or much less than the average. For example, only 
'' about, two percent of the molecules in a given volume of gas 
have speeds that exceed twice the average molecular speed. 
Thus, it should be apparent that, at high resultant speeds 
(the resultant speed being the component normal to the 
body U sin 6), the relatively small spread of molecular speeds 
predicted by the Maxwellian distribution become unimpor- 
tant in comparison with the resultant speed. The dimen- 
\ sionless number, s, the ratio of flight speed to the most prob- 
able molecular speed, which appears in free-molecule flow 
calculations, may be shown to be related to the conventional 
Mach number by a constant factor. Its significance, how- 
ever, is most logically interpreted as describing the relative 
importance of flight speed and molecular speed. 

This diminishing importance of the spread of molecular 
velocities is the reason for the constancy of the second term 
of equation (15) for values of s sin 9 greater than 2. It is 
interesting to note that the term in equation (15), which 
accounts for the translatory molecular motions, has values 
ranging from 2 kT at zero velocity to (5/2 )kT at values of s 
sin 9 of 2 and higher. This may be compared with the con- 
stant value of (3/2 )kT given in the simplified analysis of 
Sanger. The increased value at zero velocity in the present 
analysis is due to the fact that the faster-moving molecules 
are both more likely to collide with the body and also carry 
more translational energy than do the slower-speed mole- 
cules. The increase in the term with increasing values of the 
parameter s sin 9 is due, as previously inferred, to the contri- 
bution of the cross-product term of the thermal and mass 
velocities in equation (12), the unintegrated expression for the 
translational energy. The actual difference in the total 
kinetic energy term computed by means of the simplified 
theory and the more exact theory depends upon the relative 
magnitudes of the terms mU 2 /2 and i pkT. For most of the 
cases computed in this report, the term mL' 2 /2 became large 
in comparison with the term \pfcT at flight speeds above 8000 
feet per second; consequently, little difference between the 
simplified and more-exact theory would be expected at 
flight speeds greater than 8000 feet per second. 

Physically, it is clear why both the number of molecules 
striking the rear side of the body, and the translatory energy 
incident upon the rear side of the body, decrease with in- 
creasing values of s sin 6. The front of the body acts as a 
shield and prevents collisions of the rear side with all mole- 
cules except those that have absolute velocity components 
in the direction of flight with magnitude equal to the flight 
speed or higher. Consequently, as the flight speed is in- 
creased, we should expect that the incident energy on the 
rear of the plate would eventually approach zero. 




HEAT TRANSFER TO BODIES TRAVELING AT HIGH SPEED IN THE UPPER ATMOSPHERE 


•|r 

?‘ r - 4F 
• 

SKIN TEMPERATURES IN THE ABSENCE OF SOLAR RADIATION 

Figure 7 shows the magnitude of plate temperatures cal- 
culated for several altitudes in the absence of solar radiation. 



Figure 7.— The effect of flight velocity and angle of incidence on the temperature of the 
front and rear sides of an uncooled plate. Accommodation coefficient, 1.0; solar radiation 
absent; angle of incidence indicated; emissivity, 1.0; insulation between front and rear 
sides of plate. 

The marked effect of emitted radiation from the body in 
lowering the plate temperature is striking. For a large range 
of flight speeds and small plate angles, the plate tempera- 
tures are actually lower than the ambient-air temperature. 
The curves of temperature of the rear side of the body show 
an eventual decrease with flight speed in accordance with the 
discussion of the preceding paragraph. The effect upon body 
temperature of perfect thermal contact between the front 
and rear sides of a flat plate is shown in figure 8. It will be 



Figure 8.— The effect of thermal contact between front and rear sides on plate temperature. 
Accommodation coefficient, 1.0; angle of incidence indicated; emissivity, 1.0; solar radiation 
absent. 


noted that the plate temperatures are considerably lower than 
those calculated for the front side of an insulated plate, as 
shown in figure 7. It would appear that the cooling problem 
could be alleviated by insuring thermal contact between 
portions of the body inclined at positive and negative angles 
with respect to the flight path. Practically, this could be 
accomplished by boattailing the body and using a heavy 
skin of high thermal conductivity between the nose section 
and the boattailed section. An alternative scheme would 
consist of employing a circulating-type cooling system, ab- 
stracting heat from the nose section and rejecting this heat 
at the boattailed section. It is also apparent that, from a 
skin temperature standpoint, large positive body angles 
should be avoided; the body should preferably be finely 
tapered. Fortunately, this dictate coincides with aerody- 
namic requirements. 

The effect of altitude on the temperature of the front side 
of an insulated plate is also indicated in figures 7 and 8. It 
is apparent that the problem of high skin temperatures 
becomes of negligible importance above altitudes of 125 
miles for flight speeds up to 36,000 feet per second and with 
solar radiation absent. The condition of zero solar radia- 
tion corresponds, of course, to nocturnal flight. 

INTERNAL COOLING 

The effect upon body temperature of internal cooling is 
shown in figure 9 for several flight speeds and for a body 
angle of 10°. It can be seen that the body temperature 
decreases rapidly above a certain critical value of rate of 
internal cooling. This effect is due to the combined actions 
of emitted radiant energy and the rate of internal cooling. 
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At the lower rates of internal heat removal, the emitted 
radiant energy comprises the major portion of the total 
energy leaving the surface, hence changes hi the rate of 
internal cooling have little effect upon the surface tempera- 
ture. As the rate of internal cooling increases, however, 
and the surface temperature decreases, internal cooling 
becomes the controlling factor, since the emitted radiant 
energy is proportional to the fourth power of the surface 
temperature. As the rate of internal cooling is increased, a 
point is finally reached where its magnitude approaches the 
value of the incident molecular energy and when this occurs, 
the surface temperature rapidly approaches zero. 

THE EFFECT OF SOLAR RADIATION ON SKIN TEMPERATURE 

The variation of plate surface temperatures with flight 
velocity in the presence of solar radiation is shown in figure 
10. It can be seen that the effect of solar radiation on the 
surface temperatures becomes increasingly important as the 
flight altitude is increased and that the effect is strongest at 
the lower flight speeds. At altitudes greater than 100 miles, 
the presence of solar radiation markedly increases the plate 
surface temperature and this effect becomes practically inde- 
pendent of flight velocity. The implication is, of course, 
that heat transfer by molecular energy exchange is negligible. 
The strong effect of solar radiation on skin temperature at 
the higher altitudes suggests the desirability of nocturnal 
flight in order to prevent excessive skin temperatures. 

THE EFFECT OF EMISSIVITY ON SKIN TEMPERATURE 

The effect of the emissivity of the plate material on the 
surface temperatures is shown in figure 11 for the condition 
of no solar radiation and in figure 12 for the condition of max- 
imum solar radiation. It is apparent that it is desirable to 



Figure 10.— The effect of solar radiation on the surface temperature of an uncooled plate. 
Accommodation coefficient, 1.0; angle of incidence indicated; emissivity, 1.0; rear side of 
plate insulated. 
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Figure II.— Effect of emissivity on plate temperature assuming solar radiation absent. 
Accommodation coefficient, 1.0; angle of incidence indicated; rear side of plate insulated. 
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Figure 12. — Effect of emissivity on plate temperature assuming solar radiation present. 
Accommodation coefficient, 1.0; angle of incidence indicated; rear side of plate insulated. 
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utilize as high a value of skin emissivity as possible. This 
conclusion is reached because any practicable missile must 
fly over a wide range of altitude and speed and, under all 
conditions of flight speed and altitude, a high emissivity tends 
to reduce skin temperature. The effect of emissivity is most 
marked at the lower altitudes; at the higher altitudes the 
surface temperatures become almost independent of the 
emissivity. This somewhat surprising trend may be ex- 
plained as follows: At altitudes of 100 miles and higher, the 
magnitude of the molecular energy interchange is small in 
comparison with the radiant energy exchange. Thus, neg- 
lecting all energy exchange except that due to radiation, an 
energy balance on the plate may be written as «J 0 =€<rT 4 . 
The emissivity e is seen to cancel out and the plate tempera- 
ture is thus dependent only upon the magnitude of the solar 
radiation J 0 . It should be re-emphasized that these com- 
putations are valid only for gray-body radiation. 

THE EFFECT OF ACCOMMODATION COEFFICIENT ON SKIN TEMPERATURE 

The effect of the magnitude of the accommodation coef- 
ficient on the plate surface temperatures is shown in figure 
13. It may be noted that, as the accommodation coefficient 
is lowered, implying increasing difficulty of molecular energy 
exchange, the plate temperature likewise decreases. This 
phenomenon may offer a method of cooling, as it may be 
possible to decrease the accommodation coefficient by 
altering the surface characteristics of the skin. 

CONCLUSIONS 

The following conclusions can be drawn from the trends 
shown in the preceding analysis: 

1. With solar radiation absent, that is, for the condition 
of nocturnal flight, the skin cooling problem associated with 
high-speed missile flight becomes negligible at altitudes of 
125 miles and higher up to steady flight speeds of 36,000 
feet per second. 

2. The effect of solar radiation on skin temperatures 
becomes increasingly important as the flight altitude is 
increased. At an altitude of 150 miles and higher, solar 
radiation is the predominating factor that determines skin 
temperature. 

3. With solar radiation present, the effect of emissivity on 
skin temperatures becomes of decreasing importance as the 
altitude is increased. At altitudes of 125 miles and higher, 
the skin temperature at a given flight speed is practically 
independent of the emissivity and depends only on the 
intensity of the solar radiation. 

4. The strong effect of solar radiation on skin temperatures 
at altitudes above 75 miles suggests the desirability of 
nocturnal flight in order to minimize skin temperatures. 

5. In order to maintain low skin temperatures, the angle 
of inclination of the body with respect to the flight path 
should be kept as small as possible. This may be done by 
designing the body to be finely tapered and by flying at small 
angles of attack. 

6. Skin temperatures may be reduced by insuring thermal 
contact between portions of the body inclined at positive 
and negative angles with respect to the flight path. As 
much surface as possible should be inclined at negative 



Figure 13.— Variation of plate temperature with accommodation coefficient. Solar radia- 
tion absent; angle of incidence indicated; emissivity, 1.0; rear side of plate insulated. 


angles. Practically, this may be accomplished by boattailing 
the body. 

7. In order for a skin-cooling system to bo effective, the 
rate of internal cooling must be of the same order of magni- 
tude or greater than the rate at which heat is lost naturally 
by radiation from the body. Internal cooling rates below 
the natural radiation loss rate have little effect on skin tem- 
peratures. 
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8. For a practical missile designed to fly over a wide range 
of altitudes and speeds, it is desirable to make the emissivity 
of the skin as high as possible. This conclusion, however, 
is restricted to the case of skin surfaces for which the emis- 
sivity is independent of the wavelength of the emitted and 
absorbed radiant energy. 

,., 9. A decrease in accommodation coefficient results in an 
accompanying decrease in skin temperature. This phenom- 
enon may offer a method of skin cooling if it is possible to 
lower the accommodation coefficient by altering the surface 
characteristics of the skin. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., August, 1948. 
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